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We consider both massive Euler-Heisenberg-like and Euler-Heisenberg-like Electrodynamics in 
the approximation of the strong-field limit. Our analysis shows that massive Euler-Heisenberg- 
type Electrodynamics displays the vacuum birefringence phenomenon. Afterwards, we calculate 
the lowest-order modifications to the interaction energy for both classes of Electrodynamics, within 
the framework of the gauge-invariant path-dependent variables formalism. Interestingly enough, 
for massive Euler-Heisenbeg-like electrodynamics (Whichmann-Kroll) we obtain a new long-range 
(1/r®- type) correction, apart from a long-range (1/r®- type) correction to the Coulomb potential. 
However, Euler-Heisenberg-like Electrodynamics in the approximation of the strong-held limit (to 
the leading logarithmic order), display a long-range (1/r®- type) correction to the Coulomb potential. 
Again, for their non-commutative versions, the interaction energy is ultraviolet hnite. 

PACS numbers: 14.70.-e, 12.60.Cn, 13.40.Gp 


I. INTRODUCTION 

The phenomenon of vacuum polarization in Quantum Electrodynamics (QED), arising from the polarization of 
virtual electron-positron pairs and leading to nonlinear interactions between electromagnetic fields, remain as exciting 
as in the early days of QED [IHi. An example that illustrates this is the scattering of photons by photons, which 
despite remarkable progresses has not yet been confirmed Along the same line, we also recall that alternative 

scenarios such as Born-Infeld theory jl^, millicharged particles [T^ or axion-like particles [l3l - [T5l| may have more 
significant contributions to photon-photon scattering physics. 

Interestingly, it should be recalled here that the physical effect of vacuum polarization appears as a modification 
in the interaction energy between heavy charged particles. In fact, this physical effect changes both the strength and 
the structural form of the interaction energy. This clearly requires the addition of correction terms in the Maxwell 
Lagrangian to incorporate the contributions from vacuum polarization process. Two important examples of such 
a class of contributions are the Uehling and Serber correction and the Wichmann-Kroll correction, which can be 
derived from the Euler-Heisenberg Lagrangian. Incidentally, as was explained in [^ , it is of interest to notice that the 
Euler-Heisenberg result extends the Euler-Kockel calculation (in the constant background field limit), which contains 
nonlinear corrections in powers of the field strengths. Whereas the Uehling and Serber result contains corrections 
linear in the fields (but nonlinear in the space-time dependence of the background fields). We further mention that, as 
in the Euler-Heisenbeg case, Born-Infeld (BI) Electrodynamics also contains similar nonlinear corrections to Maxwell 
theory from a classical point of view, as is well known. Nevertheless, BI Electrodynamics is distinguished, since 
Bl-type effective actions arise in many different contexts in superstring theory [l^ . H^. In addition to Born-Infeld 
theory, other types of nonlinear electrodynamics have been discussed in the literature 

In this context, we also point out that extensions of the Standard Model (SM) such as Lorentz invariance violating 
scenarios and fundamental length have become the focus of intense research activity [23 - [^ . This has its origin in 
the fact that the SM does not include a quantum theory of gravitation, as to circumvent difficulties theoretical in the 
quantum gravity program. Within this context quantum field theories allowing non-commuting position operators 
has been studied by using a start product (Moyal product) (29l - l^ . Mention should be made, at this point, to a 
novel way to formulate noncommutative quantum field theory (or quantum field theory in the presence of a minimal 
length) which is implemented through a new multiplication rule which is known as Voros star-product [s^ . 

It was later shown that the physics is independent from the choice of the type of product [s^. More recently, this 
new approach has been successfully extended to black holes physics [i^ . 

Inspired by these observations, the purpose of this paper is to extend our previous studies [Hill on nonlinear 
electrodynamics to the case when vacuum polarization corrections are taken into account. The preceding studies 
were done using the gauge-invariant but path-dependent variables formalism, where the interaction potential energy 
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between two static charges is determined by the geometrical condition of gauge invariance. One important advantage of 
this approach is that it provides a physically-based alternative to the usual Wilson loop approach. Accordingly, we shall 
work out the static potential for electrodynamics which include, apart from the Maxwell Lagrangian, additional terms 
corresponding to the Uehling, massive Euler-Heisenberg-like, Euler-Heisenberg Electrodynamics in the approximation 
of the strong-held limit (to the leading logarithmic order) and for their non-commutative versions. Our results show a 
long-range ^l]j>-iyY>e correction to the Coulomb potential for both massive Euler-Heisenberg-like and Euler-Heisenberg 
Electrodynamics in the approximation of the strong-held limit (to the leading logarithmic order). Interestingly enough, 
for massive Euler-Heisenbeg-like electrodynamics (Whichmann-Kroll) we obtain a new long-range correction to 
the interaction energy. Nevertheless, for their non-commutative versions, the static potential becomes ultraviolet 
hnite. 

The organization of the paper is as follows: In Section H, we reexamine Uehling Electrodynamics in order to 
establish a framework for the computation of the static potential. In Section HI we consider Euler-Heisenberg-like 
(with a mass term) and show that it yields birefringence, compute the interaction energy for a fermion-antifermion 
pair and its version in the presence of a minimal length. In Section HI, we repeat our analysis for Euler-Heisenberg 
Electrodynamics in the approximation of the strong-held limit. Finally, in Section IV, we cast our Final Remarks. 

In our conventions the signature of the metric is (-1-1, —1, —1, —!)• 


II. BRIEF REVIEW ON THE UEHLING POTENTIAL 


As already expressed, we now reexamine the interaction energy for Maxwell theory with an additional term cor¬ 
responding to the Uehling correction (Uehling Electrodynamics). This would not only provide the setup theoretical 
for our subsequent work, but also hx the notation. To do that we will calculate the expectation value of the energy 
operator H in the physical state |$), which we will denote by We start oh our analysis by considering the 

ehective Lagrangian density [4l| : 

c = ( 1 ) 

where 


Ai (m, x) 



dr— 

T 


1 

(r-k A) 




( 2 ) 


with A = dfid^. It should be noted that Ai contains the ehect of vacuum polarization to first order in the fine 
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structure constant, a = |^, and, m, is the electron mass. 

To obtain the corresponding Hamiltonian, we will carry out the quantization of the theory using the Hamiltonian 
formalism. In this way, the canonical momenta are found to be H^ = (l — ^ AAl) and one immediately identihes 
the primary constraint n° = 0, whereas the momenta are H® = (l — -^AAi) E® (with = F^q). Accordingly, the 
canonical Hamiltonian He is 


He 


d^3 


ffa.Ao - -n.(l - ^AM) + -F, (1 - ^AAi) 


(3) 


As usual, requiring the primary constraint H*^ to be stationary, leads to the secondary constraint Ti = diJF = 0. 
It is straightforward to check that there are no further constraints in the theory. Hence, the extended Hamiltonian 
that generates translations in time then reads H = He + J d^x (uo(x)Ilo(x) -I-ui(a:)ri(a;)), where Uoix) and ui(x) 
are arbitrary Lagrange multipliers. With the aid of equation ([3]) we find that Aq (x) = [Aq (x) , H] = uq (x), which is 
an arbitrary function. Since H*' = 0 always, neither nor n° are of interest in describing the system and may be 
discarded from the theory. Therefore, the Hamiltonian becomes 


H = 


d^x 


w{x)diir 


in,(l - tAM)-'n* + (1 - ^AM) 


(4) 


where w{x) = ui(x) — Aq(x). 

According to the usual procedure, we impose a gauge condition such that the full set of constraints becomes second 
class. A convenient choice is 


i 

r2(x)= J dz''Ai,{z) = J d\x^Ai{Xx)=0. 


(5) 
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where A (0 < A < 1) is the parameter describing the spacelike straight path z® = ^* + A (a; — and ^ is a fixed point 
(reference point). There is no essential loss of generality if we restrict our considerations to = 0. The Dirac brackets 
can now be determined and we simply note the only nontrivial Dirac bracket involving the canonical variables, that 
is, 


{Ai {x ), (j/)}* = S{{x -y) - di j d\x"5^^'> (Ax - y) . (6) 

0 

In passing we also recall that the transition to a quantum theory is made by the replacement of the Dirac brackets 
by the operator commutation relations according to {A,B}* —>• (—Vr) [Al,i3]. 

With the foregoing information, we can now proceed to obtain the interaction energy. As already mentioned, in 
order to accomplish this purpose, we will calculate the expectation value of the energy operator H in the physical 
state 1$), where the physical states 1$) are gauge-invariant ones. The physical state can be written as 

1^) = |'i'(y)^(y')) = '0(y) J dz^Ai(z))-tp(y') |0), (7) 

where |0) is the physical vacuum state and the line integral appearing in the above expression is along a space-like 
path starting at y/ and ending at y, on a fixed time slice. 

Making use of the above Hamiltonian structure M , we find that 

A^(x)\'^>(y)'^|(y'))='^>(y)'^f{y')Ui{x)\0) + q^ dz^S^^'> (z - x) \<^). ( 8 ) 

Jy 

With the aid of equations (|S]) and (|1]), the lowest-order modification in a of the interaction energy takes the form 


where {H)q = 


and 


{H)^ = {H}^ + Vi + V2, 
I |0). The Vi, V 2 terms are given by 


Hi = -^ y d^x J dz^5^'^^ (x — z) J dz'^5^^^ (x — z'), 




(9) 


( 10 ) 


(11) 


where p (t) = (l + \Jl - 

We note that the term (flOl) may look peculiar, but it is nothing but the familiar Coulomb interaction plus a self¬ 
energy term [d^. Now making use of the Green function, G (z, z') = ^ ’ 1^® term m can be rewritten in 

the form 


G 2 = 


9 

a 


r dr^ r dz(vi, r dz^ciz ^)= 

JAm'^ ^ Jy Jy ^ JAm^ 


p(t) 

r |y-y'l 


( 12 ) 


Since the second and third term on the right-hand side of Eq. @ are clearly dependent on the distance between 
the external static fields, the potential for two opposite charges located at y and y/ reads 


V = - 


47r L 





(13) 


where L = |y — y^|. Accordingly, one recovers the known Uehling potential, which finds here an entirely different 
derivation. 

Before we proceed further, we wish to show that this result can be written alternatively in a more explicit form. 
Making use of [i^ 



dt 


e 


— tz 






Xn (z) 


(14) 
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we then get 


\ f 2a , , 

By the transformation, t = cosh it [i^, the functions Xn can be reduced to the form (4^ : 

Xn (z) = Kin -1 (z) - ^Kin +1 {z) - ^Ki^+s {z ), 


(15) 


(16) 


where the functions Ki denote Bessel function integrals. Hence we see that the interaction energy (with m = 1) 
becomes 


V = - 


dir L 


1 + 


2a 

Stt 


1 + 


/5L 2L^\ 


3 7 




V 2 


dtKo (t) 


l2L 


(17) 


where Ko{z) and Ki{z) are modified Bessel functions. Finally, with the aid of asymptotic forms for Bessel functions, 
it is a simple matter to find expressions for V for large and small L. 

Before concluding this subsection we discuss an alternative way of sta ting our previous result 0, which displays 
certain distinctive features of our methodology. We start by considering [421 l45l| 


V = q{Ao (0)-A(L)), 

where the physical scalar potential is given by 

Ao{t,r) = f dXAEi{t,X^)- 
Jo 

This follows from the vector gauge-invariant field expression 

Afj,{x) = (x) J dz^A^ {z) 


(18) 


(19) 


( 20 ) 


where the line integral is along a space-like path from the point ^ to x, on a fixed slice time. It is also important to 
observe that the gauge-invariant variables ([T^ commute with the sole first constraint (Gauss law), showing in this 
way that these fields are physical variables. Inasmuch as we are interested in estimating the lowest-order correction 
to the Coulomb energy, we will retain only the leading term in expression E'' = (l — ^AA4) H*. Making use of 

this last expression, equation (HU gives 


Ao{t,r) = / dXAdii- 


jo 


a 


dr 


^ J Xr JAm'^ ^ 




' dXAdJ-- 


V2-- 


( 21 ) 


Xr 


to get the last line we used Gauss law for the present theory, that is, dJE = J° (where we have included the external 
current J° to represent the presence of two opposite charges). Accordingly, for J°(t,r) = q5^^^ (r), the potential for 
a pair of static point-like opposite charges located at 0 and L, is given by 




47r L 


Stt 


nOO 

/ dr-—'—e 

Um^ T 




( 22 ) 


after substracting a self-energy term. 


III. EULER-HEISENBERG-LIKE MODEL 


Proceeding in the same way as we did in the foregoing section, we shall now consider the interaction energy for 
Euler-Heisenberg-like electrodynamics. Nevertheless, in order to put our discussion into context it is useful to describe 
very briefly the model under consideration. In such a case the Lagrangian density reads: 



1 1 


(23) 
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where have included two parameters j3 and 7 . As usual, F = Q = F^i, = — 9j/A^ and 

pnv _ j^g|. g^jgp jjiejition here that in our previous paper [T^ we have studied the domain 0 < p < 1. 

Moreover, follows from (1^51) that whenp = 2 the model contains, to order O (V/ 3 ^) and O a Euler-Heisenberg-like 

model with the appropriate identifications of the constants. Interestingly, we also observe that in the limit 7 —> 00 we 
obtain a Whichmann-Kroll model. This remark opens up the way to discuss the effect of these nonlinear corrections 
on the interaction energy, as we are going to study below. In fact, we shall consider a massive Whichmann-Kroll 
system. The motivation for this study comes from recent considerations in the context of dualities , where massive 
Born-Infeld systems play an important role. 

Having made these observations we can write immediately the field equations for p = 2: 


while the Bianchi identities are given by 


where 


r ( 1 = 0 , 


d^Ff^’' = 0 , 


F 

r = 1 -I_-— 

^2 ■ 


Also, it is straightforward to see that Gauss law becomes, 

V-D = 0, 

where D is given by 


D 


(E2-B2) (E-B) 


E-b4(E-B)B) . 


(24) 

(25) 

(26) 

(27) 

(28) 


Again, from equation (EH), for J°(t, r) = (r), we find D = ^f, where Q = This then implies that for a 

point-like charge, e, at the origin, the expression 


Q 


= I- 


E 2 


|E|, 


(29) 


tells us that, for r — 0 , the electrostatic field becomes singular at r = 0 , in contrast to the 0 < p < 1 case where the 
electrostatic field is finite. Even so, in this theory the phenomenon of birefringence is present. 

To illustrate this important feature we introduce the vectors D = dCjdFi and H = —9£/i9B: 


and 


D =r E-t-2 


B(E B) 


H = r B-2 


E(E- B) 


(30) 


(31) 


where T = 1 -f ^ (B^ — E^) — (E • B)^. We thus obtain the equations of motion 


and 


aD 

V-D = 0, —-VxH = 0, 

at 


V-B = 0, —-)-VxE = 0. 

at 


(32) 


(33) 


With the aid from equations (15(11) and (1511) we find the electric permitivity, , and the inverse magnetic permeability, 
(p“^) tensors of the vacuum, that is. 


Eij — r 



2BiBj \ 

^2 j > 


(>■-■)« =r(i., 


2EiEj\ 

—) 


(34) 
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with Di = Eij Ej and Bi = fiij Hj. 

In accordance with our previous procedure [l^, [l^, we can now linearize the above equations. To do this, it is 
advantageous to introduce a weak electromagnetic wave (Ep,Bp) propagating in the presence of a strong constant 
external held (Eo,Bo). On these assumptions, we readily hnd that, for the case of a purely magnetic held (Eq = 0), 
the vectors D and H become 


D 


1 + Sl 

2 P 


E„ 


(Ep • Bq) Be 


and 


(35) 


H = 




(Bp-Bo)Bo , 


(36) 


where we have keep only linear terms in Ep, Bp. As before, we consider the z axis as the direction of the external 
magnetic held (Bq = BqBs) and assuming that the light wave moves along the x axis, the decomposition into a plane 
wave for the helds Ep and Bp can be written as 


Ep (x, t) = , Bp (x, t) = . 


(37) 


In this case, it clearly follows that 


(^-£ 221 . 33 ) 


E 2 — 0 , 


(38) 


and 


/ 

—y — £33/^22 
\w‘‘ 


]e3=0. 


(39) 


As a consequence, we have two different situations: First, if E _L Bq (perpendicular polarization), from (15^)) E 3 = 0, 
and from (IM)l we get ^ = £ 22 /^ 33 - This then means that the dispersion relation of the photon takes the form 


n_L = 


\ 


1 + 




1 + 


3B;i 

2/32 


(40) 


Second, if E || Bq (parallel polarization), from (1551) E 2 = 0, and from (1551) we get = 633 ^ 22 - This leads to 


nil = -i/l + 


2B^ 


(41) 


Thus we verify that in the case of a generalized Euler-Heisenberg electrodynamics the phenomenon of birefringence 
is present. 

We now pass to the calculation of the interaction energy between static point-like sources for a massive Whichman- 
Kroll-like model, our analysis follows closely that of references [mil- The corresponding theory is governed by the 
Lagrangian density 


£ = — i 

4 




(42) 


Next, in order to handle the second term on the right hand in (I42L we introduce an auxiliary field ^ such that its 
equation of motion gives back the original theory. This allows us to write the Lagrangian density as 






32^2 1 ^’' 

With the redefinition 77 = 1 — equation (H51) becomes 


F^.E^’' - 


1 ..2 "m . .a 




128/32 ' 2 


(43) 


£ = - i(l - + ^A,A^. 


(44) 
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Before we proceed to work out explicitly the interaction energy, we shall first restore the gauge invariance in equation 
CT . Following an earlier procedure, we readily verify that the canonical momenta read 11^ = which results 

in the usual primary constraint 11*^ = 0, and 11* = In this way one obtains 

Hc = J d?x + i(l - 7?)"|. (45) 


The consistency condition, B = 0, leads to the constraint F = 9^11* + rn?A^ = 0. As a result, both constraints are 
second class. To convert the second class system into first class we shall adopt the procedure described previously. 
Thus, we enlarge the original phase space by introducing a canonical pair of fields 9 and Be. It follows, therefore, that 
a new set of first class constraints can be defined in this extended space: Ai = Bq + rn^O — 0 and A 2 = F + Bg = 0. 
This then shows that the new constraints are first class and, therefore, restore the gauge symmetry. From this, the 
new effective Lagrangian density, after integrating out the 9 fields, becomes 


Now, writting cr = ry + ^, the expression (HSl) can be brought to the form 

1 k / 

1 — (7 • 


1 k 

r — _ F (tF^^ _ 

^ 128 


m 

~A 


(46) 


(47) 


where k = 64/3^. 

We are now ready to compute the interaction energy. In this case, the canonical momenta are B^ = —crF°^, with 
the usual primary constraint B° = 0, and B* = crF*°. Bence the canonical Bamiltonian is expressed as 

He = / jn-a/l, + Tn» + H- A (1 I . (48) 

Time conservation of the primary constraint B° yields the secondary constraint Fi = dAF = 0. Similarly for the Va 
constraint yields no further constraints and just determines the field tr. In this case, at leading order in /3, the field a 
is given by 


cr = 




1 


3 

w 


1 




n 


2 


(49) 


which will be used to eliminate a. As before, the corresponding total (first-class) Bamiltonian that generates the time 
evolution of the dynamical variables is FI = He + f d^x (uo(x)Ilo(x) -I-ui(a:)ri(a;)), where Uo(x) and ui{x) are the 
Lagrange multiplier utilized to implement the constraints. 

In the same way as was done in the previous subsection, the expectation value of the energy operator H in the 
physical state 1$) becomes 


{H)^ = ($1 


d^x 





I5m^ 

2/32 


A 




(50) 


in this last line we have considered only quadratic terms in m^. 

In such a case, by employing dsni), the lowest-order modification in /3^ and rrA of the interaction energy takes the 
form 


{H)^ — {H)q + Vi + V2 + V3, 

where (B)g = (0| H |0). The Vi, V 2 and V 3 terms are given by 



(51) 


( 52 ) 






(53) 


V 2 = — j (fx J dzi6^^^ {x — z) J (x — z') J (x — u) J dvkS^^\:x. — v) , 

and 

1 c ;™2 4 (• py py py py 

V3 = —— J d^x J dziS^^^ (x — z) J dz'’'S^^'^ (x — z') J du^S^^'^ (x — u) J dvkS^^^ (x — v). (54) 

Finally, with the aid of the expressions (15^ . (1551) and (IMl) . the potential for a pair of static point-like opposite 
charges located at 0 and L, is given by 


V = - 




3 1 


47r L 1671/32 \ 8 Tr 


— 5m^ 


1 

Z3' 


(55) 


observe that when m = 0, the profile (1551) reduces to the known Whichmann-Kroll interaction energy. On the other 
hand, for m 0, it should be noted the key role played by the mass term in transforming the Coulomb potential 
into the Yukawa one. Interestingly enough, an unexpected feature is found. In fact, the profile (1551) displays a new 
long-range [3 correction, where its strength is proportional to m^. It is also important to observe that a analogous 
correction has been found in Born-Infeld electrodynamics in the context of very special relativity (F7j| . In this way we 
establish a new connection between nonlinear effectives theories. 

Before we proceed further, we should comment on our result. In the case of QED (Euler-Heisenberg Lagrangian 
density) the parameter 1/^2 is given by ^ = || ^ 4 ^ 7 , where me is the electron mass. In this context, we also recall 

the currently accepted upper limit for the photon mass, that is, m^ ~ 2 x 10“^® eV. Thus, for the QED case, from 
equation (1551) it follows that the second term on the right-hand side it would be detectable in long-range distances 
(^ 10® m). In other words, we see that detectable corrections induced by vacuum polarization with a mass term 
would be present at low energy scales. 

From equation (I55L it clearly follows that the interaction energy between heavy charged charges, at leading order 
in /3, is not finite at the origin. Motivated by this, one may consider the above calculation in a non-commutative 
geometry, based in findings of our previous studies Bin.' In such a case, the electric field at leading order in /3® 
and TO®, takes the form 


= 



V® 


(56) 


where it may be recalled that we are now replacing the source (x — y) by the smeared source (x — y), 

with 0 the parameter non-commutative. Now, making use of equation (HU), we readily find that 


Ao {t, r) = (t, r) -h {t, r) -H {t, r). 

The term (t,r) was first calculated in [i^, we can, therefore, write only the result: 


(i,r) = g® 


171-^9 


47 r r 


1 f°° 1 _ m^r -2 

' 2 du—^e““ 


^ , 2 du- 1 -^ 

^ ke Vu 


^47r • 


While the terms (t, r) and A)^’ (t, r), after some manipulation, can be brought to the form 

iv r.. 1 

rn du —7 


1 ( 2 ) 




du^ — P.3(3 id A 
y&' \2^ i9 ) ^ 


(t, r) = 


2^3 


3m^q 

- 

/327r®/2 


V2’4ej M2'V2’4e^ '\^2’4e^ 


(57) 

(58) 

(59) 

(60) 


where 7 ( 2 / 2 , is the lower incomplete Gamma function defined by 7 (“/&, x) = 

Inserting these expressions in equation (jl 8 p , we finally obtain the static potential for two opposite charges q located 
at 0 and L as 


V = - 


q e 


^ — mL 


47r L 
im^q^ ..j 

Jo 


'r 

VTT JT 1 ^^ y/u 


12 q^ 

/327r^/2 


y 6 ' \ 2 ^ 49 J 


du*—7® l — ^ - 7fi —^ 

u4 ^ V2>4ej y2’\2^49) ’ \2^ 49) 


(61) 


which is finite for L —>■ 0. It is a simple matter to verify that in the limit 0 —>■ 0 we recover our above result. 
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IV. LOGARITHMIC CORRECTION 


We now want to extend what we have done to Euler-Heisenberg-like electrodynamics at strong fields. As already 
mentioned, such theories show a power behavior that is typical for critical phenomena [d^. In such a case the 
Lagrangian density reads: 


C — -^F -—F log f 

^ g I 4;)^2 


(62) 


where bo and A are constants. In fact, by choosing bo = and A = , we recover the Euler-Heisenberg 

electrodynamics at strong fields [d^ . 

In the same way as was done in the previous section one can introduce an auxiliary field, to handle the logarithm 
in (1^. This leads to 


boi 


where ai = 1 —^ (1 + logC) and 02 = 32^,2 • 

A similar procedure can be used to manipulate the quadratic term in 
auxiliary field, rj, one easily finds 


£ = («! + 4a2?7) + ^« 2 - 


!L 

4 


By setting, a = ai + 4a2'rj, we then have 


1 


I 


£ = + ^j—(a - ai)A 


(63) 

Accordingly, by introducing a second 

(64) 

(65) 


6 da 2 

It is once again straightforward to apply the gauge-invariant formalism discussed in the foregoing section. The 

dL 
da 


canonical momenta read H^ = and at once we recognize the two primary constraints n° = 0 and Va = ^ = 

The canonical Hamiltonian corresponding to (1651) is 




1 / n 2 

(a - ai) 


( 66 ) 


Requiring the primary constraint n° to be preserved in time, one obtains the secondary constraint Ei = diW = 0. 
In the same way, for the constraint Fo-, we get the auxiliary field a as 


-=(l-|(l + ln?) + ^C 


1 -f 


36oB2 




2 A 2 


(i-^(i + ine) + ^e) _ 


Hence we obtain 


He = I d^x{ ^ (1 + logo 


(67) 


( 68 ) 


As before, requiring the primary constraint to be preserved in time, one obtains the auxiliary field 0 In this 

A 


case ^ = g^. Consequently, we get 

Hc= J dPx |n,a*Ao + ^ (1 + ba) . 

Eollowing the same steps that led to equation (150)1 we find that 

(F)W = ($1 Jd^x 1 $). 


(69) 


(70) 


It should be noted that this expression is similar to equation (1501) in the limit to —> 0, except by the chaged sign in 
frot of the 11"^ term. Hence we see that the potential for two opposite charges in 0 and L is given by 


r = -£i- 


1 


Ett L 6040/3^7r2 L® 


(71) 
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V. FINAL REMARKS 

Finally, within the gauge-invariant but path-dependent variables formalism, we have considered the confinement 
versus screening issue for both massive Euler-Heisenberg-like and Euler-Heisenberg Electrodynamics in the approxima¬ 
tion of the strong-field limit. Once again, a correct identification of physical degrees of freedom has been fundamental 
for understanding the physics hidden in gauge theories. Interestingly enough, their non-commutative version displays 
an ultraviolet finite static potential. The above analysis reveals the key role played by the new quantum of length in 
our analysis. In a general perspective, the benefit of considering the present approach is to provide unifications among 
different models, as well as exploiting the equivalence in explicit calculations, as we have illustrated in the course of 
this work. 
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By using the gauge-invariant, but path-dependent, variables formalism we study both massive 
Euler-Heisenberg-like and Euler-Heisenberg-like Electrodynamics in the approximation of the strong- 
held limit. It is shown that massive Euler-Heisenberg-type Electrodynamics displays the vacuum 
birefringence phenomenon. Subsequently, we calculate the lowest-order modihcations to the interac¬ 
tion energy for both classes of Electrodynamics. As a result, for the case of massive Euler-Heisenbeg- 
like electrodynamics (Whichmann-Kroll) unexpected features are found. We obtain a new long-range 
(1/r®- type) correction, apart from a long-range (1/r®- type) correction to the Coulomb potential. 
Furthermore, Euler-Heisenberg-like Electrodynamics in the approximation of the strong-held limit 
(to the leading logarithmic order), displays a long-range (1/r®- type) correction to the Coulomb 
potential. Besides, for their non-commutative versions, the interaction energy is ultraviolet hnite. 

PACS numbers: 14.70.-e, 12.60.Cn, 13.40.Gp 


I. INTRODUCTION 

The phenomenon of vacuum polarization in Quantum Electrodynamics (QED), arising from the polarization of 
virtual electron-positron pairs and leading to nonlinear interactions between electromagnetic fields, remains as exciting 
as in the early days of QED [IHi. An example that illustrates this is the scattering of photons by photons, which 
despite remarkable progress has not yet been confirmed Along the same line, we also recall that alternative 

scenarios such as Born-Infeld theory [l^, millicharged particles [T^ or axion-like particles [l3l - [T5l| may have more 
significant contributions to photon-photon scattering physics. 

Interestingly, it should be recalled here that the physical effect of vacuum polarization appears as a modification 
in the interaction energy between heavy charged particles. In fact, this physical effect changes both the strength and 
the structural form of the interaction energy. This clearly requires the addition of correction terms in the Maxwell 
Lagrangian to incorporate the contributions from vacuum polarization process. Two important examples of such 
a class of contributions are the Uehling and Serber correction and the Wichmann-Kroll correction, which can be 
derived from the Euler-Heisenberg Lagrangian. Incidentally, as explained in Q, it is of interest to notice that the 
Euler-Heisenberg result extends the Euler-Kockel calculation (in the constant background field limit), which contains 
nonlinear corrections in powers of the field strengths. Whereas the Uehling and Serber result contains corrections 
linear in the fields (but nonlinear in the space-time dependence of the background fields). We further mention that, as 
in the Euler-Heisenbeg case, Born-Infeld (BI) Electrodynamics also contains similar nonlinear corrections to Maxwell 
theory from a classical point of view, as it is well-known. Nevertheless, BI Electrodynamics is distinguished, since 
Bl-type effective actions arise in many different contexts in superstring theory [l^ . H^. In addition to Born-Infeld 
theory, other types of nonlinear electrodynamics have been discussed in the literature lii®. 

In this perspective, we also point out that extensions of the Standard Model (SM) such as Lorentz invariance 
violating scenarios and fundamental length have become the focus of intense research activity [23 - [^ . This has 
its origin in the fact that the SM does not include a quantum theory of gravitation, as to circumvent difficulties 
theoretical in the quantum gravity program. Within this context quantum field theories allowing non-commuting 
position operators has been studied by using a star product (Moyal product) j32| - [3^ . In this connection it becomes 
of interest, in particular, to recall that a novel way to formulate non-commutative quantum field theory has been 
proposed i n I38l - l40l |. The key ingredient of this development is to introduce coherent states of the quantum position 
operators |4ll | , where a modified form of heat kernel asymptotic expansion which does not suffer from short distance 
divergences has been obtained. We also point out that an alternative derivation of the coherent state approach has 
been implemented through a new multiplication rule which is known as Voros star-product [d^. Anyhow, physics 
turns out to be independent from the choice of the type of product [i^. It is worthy noting here that this type of 
non-commutativity (coherent state approach) leads to a smearing effect which is equivalent to that encountered in 
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a class of non-local theory. In other words, non-commutativity is just a sub-class of possible non-local deformation 
[3, |3| • More recently, this new approach has been successfully extended to black holes physics [3l , also in connection 
to holographic superconductors via AdS-CFT (3 |- 

Inspired by these observations, the purpose of this paper is to extend our previous studies [3 13 oa nonlinear 
electrodynamics to the case when vacuum polarization corrections are taken into account. The preceding studies 
were done using the gauge-invariant but path-dependent variables formalism, where the interaction potential energy 
between two static charges is determined by the geometrical condition of gauge invariance. One important advantage of 
this approach is that it provides a physically-based alternative to the usual Wilson loop approach. Accordingly, we shall 
work out the static potential for electrodynamics which include, apart from the Maxwell Lagrangian, additional terms 
corresponding to the Uehling, massive Euler-Heisenberg-like, Euler-Heisenberg Electrodynamics in the approximation 
of the strong-held limit (to the leading logarithmic order) and for their non-commutative versions. Our results show a 
long-range correction to the Coulomb potential for both massive Euler-Heisenberg-like and Euler-Heisenberg 

Electrodynamics in the approximation of the strong-held limit (to the leading logarithmic order). Interestingly enough, 
for massive Euler-Heisenbeg-like electrodynamics (Whichmann-Kroll) we obtain a new long-range correction to 
the interaction energy. Nevertheless, for their non-commutative versions, the static potential becomes ultraviolet 
hnite. 

The organization of the paper is as follows: In Section H, we reexamine Uehling Electrodynamics in order to 
establish a framework for the computation of the static potential. In Section HI we consider Euler-Heisenberg-like 
(with a mass term) and show that it yields birefringence, compute the interaction energy for a fermion-antifermion 
pair and its version in the presence of a minimal length. In Section HI, we repeat our analysis for Euler-Heisenberg 
Electrodynamics in the approximation of the strong-held limit. Finally, in Section IV, we cast our Final Remarks. 

In our conventions the signature of the metric is (-|-1, —1, —1, —1). 


II. BRIEF REVIEW ON THE UEHLING POTENTIAL 


As already expressed, we now reexamine the interaction energy for Maxwell theory with an additional term cor¬ 
responding to the Uehling correction (Uehling Electrodynamics). This would not only provide the setup theoretical 
for our subsequent work, but also hx the notation. To do that we will calculate the expectation value of the energy 
operator H in the physical state |$), which we will denote by (i/)$. We start oh our analysis by considering the 
ehective Lagrangian density [31 : 

c = (l - (1) 

where 


M (m, x) 



dr— 

T 


1 

(t + A) 




( 2 ) 


with A = dfj^d^. It should be noted that Ai contains the ehect of vacuum polarization to first order in the fine 
structure constant, ct = j-, and, m, is the electron mass. We may parenthetically note here that the presence of A in 
equation ([2|) does not offer problems. On the one hand, as we shall explain below we restrict to the static case. On 
the other hand, in order to compute the interaction energy, we shall make a series expansion at leading order in a. 


Before going on, two remarks are pertinent at this point. First, the modification of Coulomb’s law in equation (HD 
follows from the weak-field limit of the one-loop ehective action of quantum electrodynamics (QED). Indeed, as was 
explained in [3 |- this modihcation can be written as 


^wf = ^ J A^{x)n^u{x,y)A'' 


(3) 


where wf denotes weak held, and H^i, is the usual order-e^ polarization tensor of QED. As is well known, in momentum 
space, this tensor is given by 


A^u{k) = - k^k^) H (fc^) , 

while the momentum space spectral representation of the polarization function H {k^^ reads 

,Pii) 1 


(4) 




dt- 


/ 4m2 


t k'^ + t' 


(5) 
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with p{t) = ^1 + y/i _ i2^. Next, due to the tensor structure of n^,y(fc), C^f can then be expressed in a 

gauge invariant way, that is, in terms of F^i,. As a consequence of this, equation ([3]) reduces to the modification of 
Coulomb’s law appearing in equation (HD- We mention in passing that in | 48 l | the signature of the metric is different 
from that used in this paper. 

Second, it should be noted that the theory described by equation (HD contains higher time derivatives, hence to 
construct the Hamiltonian one must use, for example, the Ostrogradsky method |49| . Accordingly, in the theory 
under consideration the velocities have to be taken as independent canonical variables. Let us also mention here that, 
in previous studies 0113, we have shown that although theories like to equation (HD contain higher derivatives, 
in the electrostatic case the canonical momentum conjugate to velocities disappears. Hence the new Legendre 
transformation to construct the Hamiltonian reduces to the standard Legendre transformation. It should, however, 
be emphasized here that the present paper is aimed at studying the static potential of the above theory, so that A 
can be replaced by — V^. Notice that, for notational convenience, we have maintained A in equations HD and (ED, 
but it should be borne in mind that this paper essentially deals with the static case. 


Now, we move on to compute the canonical Hamiltonian. For this end we perform a Hamiltonian constraint 
analysis. The canonical momenta are found to be H'^ = (l — j^AAl) It is easy to see that n° vanishes, we 
then have the usual constraint equation, which according to Dirac’s theory is written as a weak (Ri) equation, that is, 
H*^ r; 0. It may be noted that the remaining non-zero momenta must also be written as weak equations. This leads 
to H® r: (l — ^AA4) F® (with Ei = Fio). Accordingly, the canonical Hamiltonian He is 


He 




n®5,Ao--n,(l-^AM) 'h’ 


+ i^b(l-^AAl)F®^- 


( 6 ) 


which must also be written as a weak equation. Next, the primary constraint, H*^ m 0, must be satished for all 
times. An immediate consequence of this is that using the equation of motion, Z Ri [Z,He]-, we obtain the secondary 
constraint Ti = dJV ~ 0, which must also be true for all time. In passing we recall that we are considering the static 
case, hence this new constraint does not contain time derivatives. It is straightforward to check that there are no further 
constraints in the theory. Therefore, in the case under consideration, there are two constraints, which are first class. 
According to the general theory we obtain the extended Hamiltonian as an ordinary (or strong) equation by adding all 
the first-class constraints with arbitrary constraints. We thus write H — He + J d^x (uo(a;)no(a;) -I- ui(a;)ri(x)), where 
Uo{x) and ui{x) are arbitrary Lagrange multipliers. It is also important to observe that when this new Hamiltonian 
is employed, the equation of motion of a dynamic variable may be written as a strong equation. With the aid of 
equation (0 we find that Aq (x) = [Aq (a;) ,H]=uo (x), which is an arbitrary function. Since n° Ri 0 always, neither 
A° nor H*^ are of interest in describing the system and may be discarded from the theory. In fact, the term containing 
Aq is redundant, because it can be absorbed by redefining the function 'w{x). Therefore, the Hamiltonian is now given 
as 


H = d-^ 


u;(x)a,n®--n,(l-^AM) 


-'n® + iF,,(l-^AAl)F®^- 


(7) 


where w{x) = ui{x) — Aq(x). 

It must be clear from this discussion that the presence of the new arbitrary function, w{x), is undesirable since we 
have no way of giving it a meaning in a quantum theory. Hence, according to the usual procedure, we impose a gauge 
condition such that the full set of constraints becomes second class. A convenient choice is 

1 

r2 (x) = J dz'^ Ay (z) = J dXx^Ai (Xx) = 0 . ( 8 ) 

Q* 0 

where A (0 < A < 1) is the parameter describing the space-like straight path z® = ^® -I- A (a; — and is a fixed point 
(reference point). There is no essential loss of generality if we restrict our considerations to ^® = 0. The Dirac brackets 
can now be determined and we simply note the only nontrivial Dirac bracket involving the canonical variables, that 
is. 


{A, (x) ,n^ (j/)}* = 515^'^'^ {x-y)- di J dAFJ^^^ (Xx - y). 

0 


( 9 ) 
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In passing we also recall that the transition to a quantum theory is made by the replacement of the Dirac brackets 
by the operator commutation relations according to {A, B}* —>■ {—^/h) [A, B], 

With the foregoing information, we can now proceed to obtain the interaction energy. As already mentioned, in 
order to accomplish this purpose, we will calculate the expectation value of the energy operator H in the physical 
state 1$), where the physical states 1$) are gauge-invariant ones. The physical state can be written as 

1^) = |'i'(y)^(y')) = '0(y) J ^ dzM,(z))V'(y') |o), (lo) 

where |0) is the physical vacuum state and the line integral appearing in the above expression is along a space-like 
path starting at y/ and ending at y, on a fixed time slice. The point we wish to emphasize, however, is that the 
physical fermion (y)) is not the Lagrangian fermion (ip (y)), which is neither gauge-invariant nor associate with 
an electric field. In fact, the physical fermion is the Lagrangian fermion together with a cloud (or dressing) of gauge 
fields. 

Making use of the above Hamiltonian structure M , we find that 

A^{x)\'^>{y)'^|{y'))='^>{y)'^f{y')Ui{x)\0) + q^ dziS^^'> {z - x) \<^). (11) 

Jy 

With the aid of equations dm and 0, the lowest-order modification in a of the interaction energy takes the form 

(B}^ = (B}, + V, + V2, (12) 

where (H}q = (0| B |0). The Vi, V 2 terms are given by 

Hi = ^ y d^x j dz^S^^^ (x — z) J dz'^6^^^ (x — z'), (13) 


and 


H 2 




(14) 


where p (t) = (l + y^l - 

We note that the term dT5)) may look peculiar, but it is nothing but the familiar Coulomb interaction plus a self¬ 
energy term [^. Now making use of the Green function, G (z, z') = ^ term da can be rewritten in 

the form 


2 OTT JT 




1'°^ p{t) e-v^iy-y'i 

Lm^ ^ T- |y-y'l ' 


(15) 


Since the second and third term on the right-hand side of Eq. da are clearly dependent on the distance between 
the external static fields, the potential for two opposite charges located at y and y/ reads 


y = -- 


47r L 


1 -k 


a 

Stt 


pOO 

/ dr -——e 

J Arri^ ^ 




(16) 


where L = |y — y^|. Accordingly, one recovers the known Uehling potential, which finds here an entirely different 
derivation. 

Before we proceed further, we wish to show that this result can be written alternatively in a more explicit form. 
Making use of 


we then get 


Xn (z) 




V = ( 1 + (2mL) 


47r L 


Stt'' 


(17) 


(18) 
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By the transformation, t = coshrt |54| . the functions Xn can be reduced to the form 


Xn (z) = Kin -1 {z) - ^Kin +1 {z) - ^Kin+3 {z) 


(19) 


where the functions Ki denote Bessel function integrals. Hence we see that the interaction energy (with m = 1) 
becomes 


1 


V = -^-<l + — 


47r L 


2a 


Sir 


5L 2L3 


1 + - 5 - ) Kq (2L) — — H—— ) Ki (2L) + — H—— 


3L 2L3 


dtKo {t) 


I2L 


( 20 ) 


where Ko(z) and Ki{z) are modified Bessel functions. Finally, with the aid of asymptotic forms for Bessel functions, 
it is a simple matter to find expressions for V for large and small L. 

Before concluding this subsection we discuss an alternative way of stating our previous result (TTBl) . which displays 
certain distinctive features of our methodology. We start by considering [52l l55l | 


V = q{Ao (0)-A(L)), 

where the physical scalar potential is given by 

Ao{t,r) = f dXAEi{t,X^)- 
Jo 

This follows from the vector gauge-invariant field expression 


Af,{x) = Af, (x) + dfj, J dz^Af,{z) 


( 21 ) 


( 22 ) 


(23) 


where the line integral is along a space-like path from the point ^ to x, on a fixed slice time. It is also important to 
observe that the gauge-invariant variables ( 1221 ) commute with the sole first constraint (Gauss law), showing in this 
way that these fields are physical variables. Inasmuch as we are interested in estimating the lowest-order correction 
to the Coulomb energy, we will retain only the leading term in expression i?* = (l — H*. Making use of 

this last expression, equation ( 1 ^ gives 


Ao{t,r)= I dXAdt^i-^ 


J°(Ar) 


Ar 


a dr ^ ^ 

37r T ^ ^ , 






(24) 


to get the last line we used Gauss law for the present theory, that is, diW = J° (where we have included the external 
current to represent the presence of two opposite charges). Accordingly, for J°(t,r) = (r), the potential for 

a pair of static point-like opposite charges located at 0 and L, is given by 


V = - 


47r L 


1-k 


37r 


nOO 

/ dr-—'—e 

/4m2 T 




(25) 


after subtracting a self-energy term. 


III. EULER-HEISENBERG-LIKE MODEL 


Proceeding in the same way as we did in the foregoing section, we shall now consider the interaction energy for 
Euler-Heisenberg-like electrodynamics. Nevertheless, in order to put our discussion into context it is useful to describe 
very briefly the model under consideration. In such a case the Lagrangian density reads: 





1 



(26) 


where have included two parameters (3 and 7 . As usual, T = Q = F^i, = and 

p^iv _ g^jgp jjiejition here that in our previous paper [T^ we have studied the domain 0 < p < 1 . 

Moreover, follows from (1^51) that whenp = 2 the model contains, to order O (V/ 3 ^) and O (Vt^); a Euler-Heisenberg-like 
model with the appropriate identifications of the constants. Interestingly, we also observe that in the limit j —>■ 00 we 
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obtain a Whichmann-Kroll model. This remark opens up the way to discuss the effect of these nonlinear corrections 
on the interaction energy, as we are going to study below. In fact, we shall consider a massive Whichmann-Kroll 
system. The motivation for this study comes from recent considerations in the context of dualities , where massive 
Born-Infeld systems play an important role. 

Having made these observations we can write immediately the field equations for p = 2: 


p ( _ 


= 0 , 


while the Bianchi identities are given by 


= 0 , 


where 

r = iH- - —. 

Also, it is straightforward to see that Gauss law becomes, 

V-D = 0, 


where D is given by 


D = 


(E 2 _ B 2 ) 

^ 2^2 


/32q2 


E 


(E-B)B 


(27) 

(28) 

(29) 

(30) 

(31) 


Again, from equation (l30l) . for J°(<, r) = (r), we find D = ^f, where Q = This then implies that for a 

point-like charge, e, at the origin, the expression 


Q 

J .2 



(32) 


tells us that, for r —>• 0, the electrostatic field becomes singular at r = 0, in contrast to the 0 < p < 1 case where 
the electrostatic field is finite. Even so, in this theory the phenomenon of birefringence is present. Before going into 
details, we would like to recall that birefringence refers to the property that polarized light in a particular direction 
(optical axis) travels at a different velocity from that of light polarized in a direction perpendicular to this axis. 
Indeed, due to quantum fluctuations the QED vacuum has this property, as we are going to show. 

To illustrate this important feature we introduce the vectors D = 9£/9E and H = —(9£/9B; 


and 


D = r 


(^E-t-2 


B(E-B)\ 

7^ 


H = r 


^B -2 


E(E- B)^^ 


(33) 


(34) 


where T = 1 -|- ^ (B^ — E^) — (E • B)^. We thus obtain the equations of motion 


and 


V-D = 0, 


5D 


- Vx H = 0, 


(35) 


V - B = 0, 


— + V X E = 0. 
ot 


(36) 


With the aid from equations (l33l) and (l34l) we find the electric permittivity, £ij , and the inverse magnetic permeability, 
(p,“^)^^., tensors of the vacuum, that is. 


£ij — r 



2BiBj \ 
7 ^ ) ’ 


G')«=r(«, 


2EiEj\ 

—) 


(37) 
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with Di = Eij Ej and Bi = fiij Hj. 

In accordance with our previous procedure [l^, [l^, we can now linearize the above equations. To do this, it is 
advantageous to introduce a weak electromagnetic wave (Ep,Bp) propagating in the presence of a strong constant 
external held (Eo,Bo). On these assumptions, we readily hnd that, for the case of a purely magnetic held (Eq = 0), 
the vectors D and H become 


D = 1 


3. 

2/32 


E. 


P ^ 2 

T 


(Ep • Bq) Be 


and 




Bp + 




(Bp • Bq) Bq 


(38) 


(39) 


where we have keep only linear terms in Ep, Bp. As before, we consider the z axis as the direction of the external 
magnetic held (Bq = BqBs) and assuming that the light wave moves along the x axis, the decomposition into a plane 


wave for the helds Ep and Bp can be written as 


Ep (x, t) = , Bp (x, t) = . 

In this case, it clearly follows that 

/fc 2 


and 


( — £ 22/^33 ] E 2 — 0, 


[ fc2 

—y — £33/^22 
\w‘‘ 


]e3=0. 


(40) 


(41) 


(42) 


As a consequence, we have two different situations: First, if E _L Bq (perpendicular polarization), from E 3 = 0, 
and from (l4T|l we get = £ 224 * 33 - This then means that the dispersion relation of the photon takes the form 


n_L = 


\ 


1 + 




1 + 


3B;i 

2/32 


(43) 


Second, if E || Bq (parallel polarization), from (HD) E 2 = 0, and from (11^ we get = £ 334 * 22 - This leads to 


nil = -i/l + 


2 B^ 


(44) 


Thus we verify that in the case of a generalized Euler-Heisenberg electrodynamics the phenomenon of birefringence 
is present. 

We now pass to the calculation of the interaction energy between static point-like sources for a massive Whichman- 
Kroll-like model, our analysis follows closely that of references [3 El- The corresponding theory is governed by the 
Lagrangian density 


£ = _ -I_ — 

4 32^2 




(45) 


Next, in order to handle the second term on the right hand in (I45L we introduce an auxiliary field ^ such that its 
equation of motion gives back the original theory. This allows us to write the Lagrangian density as 






32^2 1^’' 

With the redefinition 77 = I — equation (H51) becomes 


F^.F^’' - 


1 ..2 "m . .a 




128/32'’ ' 2 


(46) 


£ = - i(l - + ^A,A^. 


(47) 
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Before we proceed to work out explicitly the interaction energy, we shall first restore the gauge invariance in equation 
m- Following an earlier procedure, we readily verify that the canonical momenta read 11^ = which results 

in the usual primary constraint 11*^ = 0, and 11* = In this way one obtains 

Hc = J d^x + i(l -rj^Y (48) 


The consistency condition, n = 0, leads to the constraint F = 9^11* + rn?A^ = 0. As a result, both constraints are 
second class. To convert the second class system into first class we shall adopt the procedure described previously. 
Thus, we enlarge the original phase space by introducing a canonical pair of fields 9 and Be. It follows, therefore, that 
a new set of first class constraints can be defined in this extended space: Ai = IIo + rn^O — 0 and A 2 = F + Ilg = 0. 
Notice that this new 0-field is not to be confused with the usual non-commutative parameter. This then shows that 
the new constraints are first class and, therefore, restore the gange symmetry. As is well known, this procedure 
reproduces the usual Stiickelberg formalism. From this, the new effective Lagrangian density, after integrating out 
the 6 fields, becomes 


r =_p 

// 




2 

Now, writing a = rj + ^, the expression (14^ can be brought to the form 


(49) 


(50) 


where k = 64,5^. 

We are now ready to compute the interaction energy. In this case, the canonical momenta are 11^ = —crF°^, with 
the usual primary constraint 11° = 0, and 11* = tTF*°. Hence the canonical Hamiltonian is expressed as 

iJc = y'd"a;|n*5,Ao + —+ |. (51) 

Time conservation of the primary constraint n° yields the secondary constraint Fi = dAF = 0. Similarly for the Va 
constraint yields no further constraints and just determines the field a. In this case, at leading order in /3, the field a 
is given by 


cr = 




1 


2/32 


A 2/32 




(52) 


which will be used to eliminate a. As before, the corresponding total (first-class) Hamiltonian that generates the time 
evolution of the dynamical variables is FI = Hq + / (ito(cc)no(a;) -I-ui(a:)ri(a;)), where Uo(x) and ui{x) are the 
Lagrange multiplier utilized to implement the constraints. 

In the same way as was done in the previous subsection, the expectation value of the energy operator FI in the 
physical state |$) becomes 


{ H )^ = ($1 


d^x 





I5m^ 

2/32 


A 




(53) 


in this last line we have considered only quadratic terms in w?. 

In such a case, by employing dSSl), the lowest-order modification in /3^ and m? of the interaction energy takes the 
form 


{H)^ = (H)o + + F2 + F3, 

where (iL)o = (0| H |0). The Vi, V 2 and V 3 terms are given by 



(54) 


(55) 
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^2 = - 


15g4 

8 ^ 


dh 


dziS^^'^ (x — z) / dz'^S^^^ (x — z') / du^S^^^ (x — u) / dvkS^^^ (x — v), 


(56) 


and 

1 c;™2 4 (• py py py py 

V3 = —— J d^x J dziS^^^ (x — z) J dz'’'S^^'^ (x — z') J du^S^^'^ (x — u) J dvkS^^^ (x — v). (57) 

Finally, with the aid of the expressions (1551) . (I5S1) and (IS3, the potential for a pair of static point-like opposite 
charges located at 0 and L, is given by 


V = - 




3 1 


47r L 1671/32 \8Tr 


— 5m^ 


1 

Z3' 


(58) 


observe that when m = 0, the profile (1551) reduces to the known Whichmann-Kroll interaction energy. On the other 
hand, for m 0, it should be noted the key role played by the mass term in transforming the Coulomb potential 
into the Yukawa one. Interestingly enough, an unexpected feature is found. In fact, the profile (1551) displays a new 
long-range [3 correction, where its strength is proportional to m^. It is also important to observe that a analogous 
correction has been found in Born-Infeld electrodynamics in the context of very special relativity (s^ . In this way we 
establish a new connection between nonlinear effectives theories. 

Before we proceed further, we should comment on our result. In the case of QED (Euler-Heisenberg Lagrangian 
density) the parameter 1/^2 is given by ^ = || ^ 4 ^ 7 , where me is the electron mass. In this context, we also recall 

the currently accepted upper limit for the photon mass, that is, m^ ~ 2 x 10“^® eV. Thus, for the QED case, from 
equation (1551) it follows that the second term on the right-hand side it would be detectable in long-range distances 
(^ 10® m). In other words, we see that detectable corrections induced by vacuum polarization with a mass term 
would be present at low energy scales. 

From equation (1551) . it clearly follows that the interaction energy between heavy charged charges, at leading order 
in /3, is not finite at the origin. Motivated by this, one may consider the above calculation in a non-commutative 
geometry, based in findings of our previous studies Bin.' In such a case, the electric field at leading order in /3® 
and TO®, takes the form 


= 



V® 


(59) 


where it may be recalled that we are now replacing the source (x — y) by the smeared source (x — y), 

with 0 the parameter non-commutative. Now, making use of equation dUD, we readily find that 


Ao {t, r) = (t, r) -h {t, r) -h {t, r). 

The term (t,r) was first calculated in [s^, we can, therefore, write only the result: 


(i,r) = g® 




47 r r 


1 f°° 1 _ m^r-2 

' 2 du—^e““ 


^ , 2 du- 1 -^ 

^ ke Vu 


^47r • 


While the terms (t, r) and A)^’ (t, r), after some manipulation, can be brought to the form 

IV r , 

rn du 


1 ( 2 ) 




y&' \2^ i9) ^ 


(t, r) = 


2^3 


3m^q 

- 

/327r®/2 


V2’4ej M2'V2’4e^ '\^2’4e^ 


(60) 

(61) 

(62) 

(63) 


where 7 ( 2 / 2 , is the lower incomplete Gamma function defined by 7 (“/&, x) = “. 

Inserting these expressions in equation (EH, we finally obtain the static potential for two opposite charges q located 
at 0 and L as 


V = - 


q e 


^ — mL 


47r L 
im^q^ ..j 

Jo 


'r 

VTT JT 1 ^^ y/u 


12q^ 

/327r^/2 


y6 ' \2^ 46 J 


du*—7® l — ^ - 7fi —^ 

u4 ^ V2>4ej y2’\2^4e) ’ \2^ 46) 


(64) 


which is finite for L —>■ 0. It is a simple matter to verify that in the limit 0 —>■ 0 we recover our above result. 
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IV. LOGARITHMIC CORRECTION 


We now want to extend what we have done to Euler-Heisenberg-like electrodynamics at strong fields. As already 
mentioned, such theories show a power behavior that is typical for critical phenomena [ssj] . In such a case the 
Lagrangian density reads: 


C — -^F -—F log f 

^ g I 4;)^2 


(65) 


where bo and A are constants. In fact, by choosing bo = and A = , we recover the Euler-Heisenberg 

electrodynamics at strong fields [s^. 

In the same way as was done in the previous section one can introduce an auxiliary field, to handle the logarithm 
in (p51) . This leads to 


where ai = 1 —^ (1 + logC) and 02 = 32 ;, 


boi 


( 66 ) 


A similar procedure can be used to manipulate the quadratic term in (16611 . Accordingly, by introducing a second 
auxiliary field, rj, one easily finds 

(67) 


£ = («! + 4a2?7) + ^« 2 - 


!l 

4 


By setting, a = ai + 4a2'rj, we then have 


1 


I 


£ = + ^j—(a - ai)A 


( 68 ) 


6402 

It is once again straightforward to apply the gauge-invariant formalism discussed in the foregoing section. The 

dL 
da 


canonical momenta read H^ = and at once we recognize the two primary constraints n° = 0 and Va = ^ = 

The canonical Hamiltonian corresponding to (1681) is 


Hc= I d^x 


-B^ - 
2 6402 


1 / n 2 

(cr - oi) 


(69) 


Requiring the primary constraint n° to be preserved in time, one obtains the secondary constraint Ei = diW = 0. 
In the same way, for the constraint Fo-, we get the auxiliary field a as 




1 -f 


36oB2 




2 A 2 


(i-^(i + ine) + ^e) _ 


Hence we obtain 


He = I d^x{ ^ (1 + logo 


(70) 


(71) 


As before, requiring the primary constraint to be preserved in time, one obtains the auxiliary field 0 In this 
case ^ Consequently, we get 


He = j d^x{ n,a*Ao + ^ (1 + h) 


(72) 


Eollowing the same steps that led to equation (l53)l we find that 


(F)W = ($1 J d^x 1 $). (73) 

It should be noted that this expression is similar to equation (1531) in the limit m —> 0, except by the changed sign in 
front of the 11^ term. Hence we see that the potential for two opposite charges in 0 and L is given by 


r = -£l- 


1 


47r L 6040/3^7r2 L® 


(74) 
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V. FINAL REMARKS 

Finally, within the gauge-invariant but path-dependent variables formalism, we have considered the confinement 
versus screening issue for both massive Euler-Heisenberg-like and Euler-Heisenberg Electrodynamics in the approxima¬ 
tion of the strong-field limit. Once again, a correct identification of physical degrees of freedom has been fundamental 
for understanding the physics hidden in gauge theories. Interestingly enough, their non-commutative version displays 
an ultraviolet finite static potential. The analysis above reveals the key role played by the new quantum of length 
in our analysis. In a general perspective, the benefit of considering the present approach is to provide a unification 
scenario among different models, as well as exploiting the equivalence in explicit calculations, as we have illustrated 
in the course of this work. 
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